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1. Introduction 


Gallé, R., Vesztergom, N. & Somogyi, T. 2011: Environmental conditions affect- 
ing spiders in grasslands at the lower reach of the River Tisza in Hungary. — 
Entomol. Fennica 22: 29-38. 


The aims of the present study were to reveal the differences between grassland 
types, and to identify the local and landscape parameters that influence the spider 
assemblages at the reach of the River Tisza. The rarefied number of spider species 
was negatively correlated with the proportion of forests in a radius of 500 m 
around each site. A negative correlation was found between the number of grass- 
land specialist species and the proportion of the forests, but the number of forest 
species increased significantly with the neighboring forest area. The relative area 
of neighboring forests, the number of plant species and regular flooding played 
major roles in shaping of the species composition of spiders. The results of the 
present study emphasize the importance of the effect of habitat landscape proper- 
ties on spider assemblages. The structural diversity of the landscape may influ- 
ence species richness and composition of the habitats. 
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habitats (Ward et al. 2002, Lambeets et al. 2008, 
2009). 


In the Hungarian Great Plain, grasslands are 
threatened by the expansion of agricultural activ- 
ity (Horvath et al. 2009). During the river regula- 
tions in the 19" century the historical floodplain 
was divided by dikes into floodplain and non- 
flooded, protected areas. The remnants of the pro- 
tected grasslands are mainly used as pastures or 
meadows. 

These small patches are enclosed between 
highly modified landscape elements, i.e. arable 
fields and forest plantations of mainly non-native 
tree species. The natural disturbance of the regu- 
lar flooding events and the related changes in the 
vegetation structure are expected to affect struc- 
tural and microclimatic parameters in floodplain 


The distribution of organisms across such a 
mosaic landscape is influenced by numerous fac- 
tors and presumably by the complex relationships 
among these factors (Turner 2005). 

Among the habitat features, the composition 
and structure of vegetation (Ysnel & Canard 
2000, Heikkinen & MacMahon 2004) influence 
the species composition and diversity of the in- 
vertebrate assemblages. Plant species richness 
and the height of the vegetation is supposed to 
have a major influence on the spider assemblages, 
as higher plant species richness offers more di- 
verse habitat structure and more potential sites for 
web building (Jeanneret et al. 2003a, Horvath et 
al. 2009). 
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On the landscape scale, heterogeneity of the 
surrounding habitat patches may influence the di- 
versity and the number of species of the assem- 
blages at a given patch (Duelli 1997). Although 
numerous authors have demonstrated the signifi- 
cance of the spatial structure of the landscape on 
invertebrate assemblages (e.g. Diekétter 2008, 
Ricketts 2001), it is difficult to generalize the ef- 
fect of the landscape heterogeneity, as various 
taxa react differently to the habitat and landscape 
features, since species ecology and dispersal abil- 
ities are different for every organism (Burel & 
Baudry 1995, Jeanneret et al. 2003a). 

The present study aimed at testing the effect 
of habitat parameters (soil humidity, regular 
flooding, structure and cover of the vegetation) 
and landscape (proportion of grasslands, forests 
and arable fields) on the spider assemblages, and 
to reveal the differences between assemblages of 
the floodplain and on protected grasslands. Spi- 
ders were chosen for the study as they are the 
most diverse and abundant predatory invertebrate 
group of grasslands (Wise 1993). 


The following hypotheses were tested: 


(1) In terms of the habitat features we expected 
the number of species to increase with more 
structured vegetation, as it is well-known that 
vegetation height and density is important for 
spiders (e.g. Greenstone 1984). We did not 
expect the number of species to change with 
regular flooding event as floodplain habitats 
do not necessarily harbor more species but 
rather exhibit different species compositions 
and different diversity patterns (Lambeets et 
al. 2009). 

(2) According to Entling et al. (2007) and Lam- 
beets et al. (2008, 2009) we also expected the 
vegetation structure, soil moisture and regular 
flooding to alter the species composition of 
spiders. 

(3) As regards the effect of the landscape compo- 
sition, we expected the number of species to 
increase with increasing proportion of forests 
around the study sites, as a consequence of the 
occurrence of forest specialist and generalist 
species in the grasslands (Usher et al. 1993). 
We also expected the number of species to in- 
crease with the increasing proportion of 
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grasslands around the study sites, as the the- 
ory of island biogeography states that more 
isolated habitats retain less habitat specialist 
species than less isolated ones (e.g. Lövei et 
al. 2006). 

(4) According to the above-mentioned hypothe- 
ses we expected both the proportion of forests 
and grasslands to influence the species com- 
position of spider assemblages. 


2. Materials and methods 
2.1. Study sites and sampling 


Two habitat complexes of the lower Tisza-valley 
were selected for sampling spider assemblages. 
The southern habitat complex (Vesszés) is situ- 
ated near Szeged (46°17°30”N; 20°14’20’E), 
where the landscape mainly consists of arable 
fields with small patches of grasslands, and oak 
and poplar forest plantations are embedded in the 
matrix of arable fields. The Doc habitat complex 
lies 40 km north of Szeged (46°26’20”N; 20°10’ 
40”E), in a structurally more complex landscape. 
Higher proportion of semi-natural grasslands and 
forests occur here. 

The spider assemblages were studied at six 
grasslands (sampling sites) in the southern habitat 
complex and four in the northern habitat com- 
plex, respectively. In each habitat complex, hay- 
meadow, saline grassland, floodplain meadow, 
degraded grassland, and at the southern habitat 
complex a pair of dike-side meadows (i.e. flood- 
plain and protected) were studied. 

To sample the fauna pitfall traps were applied 
(diameter 65 mm, filled with ethylene glycol as 
preservative, Koivula 2003, Schmidt et al. 2006). 
At each site three lines of traps (sampling plots) 
were placed, each line consisted of five traps 
keeping a distance of 4 m between them. The av- 
erage distance between the lines was 150 m. Prior 
to the data analysis, the data of the five traps were 
pooled. We expected an underestimation of the 
abundance of web-building species, as pitfall 
traps measure the activity-density of species at 
the ground level, thus they capture the spiders 
with an active hunter lifestyle on the soil surface 
more efficiently (Topping & Sunderland 1992, 
Urak & Samu 2008). However, pitfall traps offer 
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a relatively good alternative for comparing the 
activity-density of epigaeic spider fauna of differ- 
ent habitats (Scmidt et al. 2005, Oberg et al. 
2007). The traps were kept open for four 2-week 
sampling periods (27 May—13 June, 30 June—15 
July, 16-30 Aug., 18 Sept.—01 Oct.) in 2007 to 
cover the main activity period of the spider spe- 
cies. The data of the four sampling periods were 
pooled. 

To characterize the structure of the vegeta- 
tion, the number of plant species and the total 
cover of the vegetation at the ground level, at 10 
and 40 cm above the ground, were measured 
within 1 x 1 meter quadrates near the traps, as 
structurally complex habitats may provide more 
potential sites for web building, and because spi- 
ders are strongly influenced by the habitat struc- 
ture (Schwab et al. 2002). 

To assess the features of the landscape, the 
proportion of land-use types was measured in a 
radius of 500 m around each site, based on aerial 
photographs (i.e. grasslands, forests and arable 
fields), as several studies suggest that landscape 
composition at scales of 500 m radius can be rele- 
vant for spiders (Clough et al. 2005, Schmidt et 
al. 2005). The remaining landscape elements 
(e.g. surface of the river) were not entered into the 
analysis in order to decrease the number of fac- 
tors which were tested. 


2.2. Data analysis 


Habitat structure affects the trappability of spider 
individuals, which may result in biased data for 
studies comparing species richness of habitats 
with different structures (Melburne 1999). Due to 
the different habitat properties of the sampling 
plots and the large variation in the number of indi- 
viduals recorded over the trapping period, rar- 
efaction was used to standardise the number of 
species recorded within each of the sampling 
plots (Heck etal. 1975, Gotelli & Colwell 2001). 
As regards hypothesis | and 3, the univariate 
regression tree method (URT) was used to gain 
insight into the structure of the data and identify 
the interactions between the variables (Tree 
Package, Ripley 2009). The method performs a 
binary recursive partitioning of the dataset. It also 
offers the opportunity to identify the influential 
explanatory variables (Crawley 2007). 


To assess the significance of their effect, the 
explanatory variables suggested by the URT me- 
thod were subjected to a linear mixed-effect 
model with nested random effects using the li- 
brary nlme in R (Pinheiro et al. 2007). The effect 
of the different habitat complexes and sampling 
sites were used as random effects and the explan- 
atory variables suggested by the regression tree 
method were used as fixed effects in the mixed ef- 
fect linear model (Crawley 2007). The collected 
species were categorized according to their habi- 
tat requirements (Buchar & Ruzicka 2002, Szita 
et al. 2006, Schmidt & Hanggi 2007, Batary et al. 
2008). The significant landscape scale parame- 
ters were included in two subsequent linear mod- 
els to assess their effects on the number of forest 
and grassland specialist species. 

Rényi’s diversity ordering was used to com- 
pare the diversity of the floodplain and protected 
habitats (hypothesis 1). Diversity ordering pro- 
vides the diversity profiles of the assemblages. 
These curves consist of a series of diversity indi- 
ces, as the scale parameter is changed, including 
indices sensitive to both rare and abundant spe- 
cies (Carranza et al. 2007). An assemblage is con- 
sidered to be more diverse than another if its pro- 
file lies above and they do not intersect (Toth- 
mérész 1995). R software (R Development Core 
Team 2007) with the BiodiversityR Package 
(Kindt 2008) was used for the calculations. 

The influence of the measured explanatory 
variables was tested on spider assemblages (hy- 
pothesis 2 and 4) by using the multivariate regres- 
sion tree (MRT) method. This method is an exten- 
sion of univariate regression trees having multi- 
variate response (De’ ath 2002). MRT does not as- 
sume particular form of relationship between spe- 
cies abundances and explanatory variables. This 
method can be used to identify interactions be- 
tween explanatory variables (De’ath 2002, 
2010). 

The size of the tree was selected based on the 
minimum tree size that falls below 1 SE of the 
minimum cross-validation estimate. The final 
tree size was selected based on the frequently oc- 
curring number of leaves from 100 individual 
trees (Work et al. 2004). 

The species that characterize each leaf and 
node of the MRT were identified using the indica- 
tor species index (IndVal) based on the relative 
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Forest (%)< 37.71 Forest (%)>37.71 
| 











Veg.cov. (%)<90.83 | Veg.cov. (%)>90.83 





=9 





Arable field (%)< 32.41| Arable field (%)> 32.41 z 
n=7 


n=6 n=8 

Fig. 1. Univariate regression tree for the rarefied spe- 
cies richness. Node 1 is based on Forest (%), the pro- 
portion of forests around the sampling plots; node 2 
on Veg. cov. (%), the total vegetation cover; node 3 on 
Arable field (%), the proportion of the arable field. The 
different lengths of the lines following each split are 
proportional to the variance explained by the split. The 
numbers below each leaf show the number of samp- 
les falling into that leaf. 


abundance and relative frequency of occurrence 
(Dufrene & Legendre 1997, De’ath 2002). 


3. Results 


A total number of 3547 spider specimens were 
collected, representing 95 species in 17 families 
(Appendix, downloadable from http://ojs.tsv.fi/ 
index.php/entomolfennica/). The species rich- 
ness of grassland patches ranged between 23 and 
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34, thus the species composition differed consid- 
erably between the studied grasslands. 


3.1. Hypothesis 1: The influence of habitat 
features on the rarefied number of the 
species 


According to the univariate regression tree me- 
thod the total vegetation cover is the only habitat 
parameter that had a possible influence on the rar- 
efied number of species in case of higher propor- 
tion of forests (Fig. 1). However, the vegetation 
cover had no significant effect on the rarefied 
number of the species according to the results of 
the linear mixed-effect model (¢,, = —0.361; p = 
0.722). 

The diversity profiles of the spider assem- 
blages inhabiting the floodplain and protected 
habitats are given in Fig 2. The curve for the pro- 
tected grasslands runs above that for the 
floodplain grasslands and the two curves do not 
intersect. Thus, the assemblages of the protected 
grasslands were more diverse for the whole range 
of the scale parameter. 


3.2. Hypothesis 2. The influence of habitat 
features on the species composition 


The cross validation of 100 MRT trees yielded 
trees of two, four and five leaves, 11, 46 and 40 
times, respectively, with three trees resulting in 
6-8 leaves. The four-leaf tree is presented in Fig 
3. The MRT, indicating the influence of the num- 


— protected grasslands 


— floodplain grasslans 


Fig. 2. The results of 
Rényi’s diverisity order- 
ing. The diversity pro- 
files do not intersect, 
thus the assemblages 
of the protected grass- 
lands are more diverse. 


3 4 
Scale parameter 
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No. of plant 
species >=13 


No. of plant 
species <13 


Node 2 
T are (0.527) * 





H rad (0.749) ** 
M ros (0.555) ** 
G luc (0.547) * 
G ruf (0.555) * 
Z seg (0.444) * 
Z lon (0.514) * 


n=9 









n=12 


Forest (%)< 37.71, Forest (%)>=37.71 






Protected | Floodplain 


O api (0.994) *** 
A ele (0.666) ** 
A leo (0.666) ** 
A mac (0.825) ** 
E den (0.578) ** 
P lat (0.679) ** 
A hum (0.5) * 
n=6 D lut (0.476) * 
P ame (0.333) * 
P cri (0.585) * 

P deg (0.550) * 
P pro (0.602) * 

T str (0.512) * 

Z mun (0.5) * 












Node 3 






D con (0.637) ** 
P ala (0.878) ** 
Z lat (0.512) * 
O sim (0.5) * 


n=3 







Fig. 3. Multivariate regression tree for spiders. Bray-Curtis dissimilarity was used for splitting. The final tree size 
was selected based on the most frequently occurring number of leaves from 100 individual trees. The numbers 
below each leaf show the number of samples falling into that leaf. The abbreviated names of species with signifi- 
cant indicator value are given at each leaf, the indicator values are given in parenthesis (*** p<0.001, ** p<0.01, * 
p<0.05). Abbreviations: T are: Thanatus arenarius, H rad: Hogna radiata, M ross: Micaria rossica, G luc: Gna- 
phosa lucifuga, Z seg: Zelotes segrex, Z lon: Z. longipes, D con: Diplostyla concolor, P ala: Pardosa alacris, Z lat: 
Z. latreillei, O sim: Ozyptila simplex, O api: Oedothorax apicatus, A ele: Antistea elegans, A leo: Arctosa leopar- 
dus, E den: Erigone dentipalpis, P lat: Pirata latitans, A. hum: Araeoncus humilis, D lut. Drassyllus lutetianus, 

P ame: P. amentata, P cri: P. cribrata, P deg: Pachygnatha degeeri, P pro: Pardosa proxima, T str. T. striatus, 


Z mun: Z. mundus. 


ber ofplant species (node 2, Fig. 3), accounted for 
12.9% of the total variance and the regular flood- 
ing events (node 3, Fig. 3) accounted for 8.9% of 
the total variance. 

The IndVal analysis identified one indicator 
species of the first leaf (Thanatus arenarius Tho- 
rell, 1872). This species lives in open, dry habi- 
tats. 

For the second leaf the IndVal analysis identi- 
fied 6 species. These spiders occur in the natural 
alkali grasslands of the Great Hungarian Plain 
(Szita et al. 2006, Batary et al. 2008). 

For third leaf the IndVal analysis identified 
two forest specialist species (Pardosa alacris 
(C.L. Koch, 1833) and Diplostyla concolor (Wi- 
der, 1834)), and the generalist Zelotes latreillei 
(Simon, 1878) and Ozyptila simplex (O.P.-Cam- 
bridge, 1862) (Buchar & Ruzicka 2002). 


For the fourth leaf we identified 14 indicator 
species. Out of these spiders the generalist Oedo- 
thorax apicatus (Blackwall, 1850), Erigone den- 
tipalpis (Wider, 1834) and Pachygnatha degeeri 
Sundevall, 1830 occur under the regularly dis- 
turbed conditions of agricultural fields and 11 
species are associated with wet meadows and 
pond margins (Buchar & Ruzicka 2002, Schmidt 
& Hanggi 2007). 


3.3. Hypothesis 3. The effect of landscape 
scale parameters on the rarefied number 
of the species 


According to the URT the two most influential 
parameters were the proportion of forests and ara- 
ble fields (Fig. 1). As predicted, the proportion of 
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forests had a significant effect according to the 
subsequent GLMM (¢,, = —2.228, p = 0.039). 
However, no significant relationship was found 
between the proportion of arable fields and the 
rarefied number of the species (¢,, = —1.801, p = 
0.089). The subsequent linear models revealed 
the significant negative effect of forests on the 
number of grassland specialist species (t,, = — 
2.166, p = 0.043) and the marginally significant 
positive effect on the forest specialist species (¢,, 
= 2.024, p = 0.057). 


3.4. Hypothesis 4. The influence of landscape 
scale parameters on the composition 
of species 


The MRT analysis identified only one influential 
landscape scale parameter. The proportion of fo- 
rests in a radius of 500 m around each site (Fig. 3, 
node 1) accounted for 26.2% of the total variance. 


4. Discussion 


4.1. No influential habitat parameter on the 
rarefied number of the species, but there is 
a higher diversity at the protected grasslands 


Although numerous studies have demonstrated 
the effect of plant species richness and habitat 
structure on the number of species and the diver- 
sity of spiders (e.g. Jeanneret et al. 2003b, Tews et 
al. 2004), we found no significant relationship 
between the vegetational parameters and the rar- 
efied number of the species. Gallé & Torma 
(2009) and Gallé et al. (2010) demonstrated that 
different habitat types may exhibit spider assem- 
blages with similar number of species but differ- 
ent species compositions. Thus, for a detailed 
study into the relationship between the spiders 
and the habitat parameters the examination of the 
species abundance data is required in order to as- 
sess the impact of these parameters on the spider 
assemblages without the loss of information 
when summarizing the data on the assemblages in 
a single value such as species richness (Jeanneret 
et al. 2003b, c). 

In accordance with Bell et al. (1999) and 
Bonn et al. (2002), our results also indicate that 
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water regime play an important role in the shap- 
ing of the diversity of spider assemblages at rive- 
rine landscapes. Numerous studies have demon- 
strated that spiders are particularly sensitive to 
flooding. The spider diversity is not only influ- 
enced by the occurrence of flooding events but 
also by the duration and frequency of these events 
(Uetz 1976, Paetzold et al. 2008). 


4.2. The number of plant species 
and the regular flooding events influenced 
the species composition of spiders 


The vegetation and flooding events are con- 
straints for spider occurrence, indicating the im- 
portance of the quality of local habitats. In low- 
land floodplains, flooding dynamics determine 
the composition of the species and the structure of 
the vegetation which in turn influence the compo- 
sition of spider fauna (Ballinger et al. 2005). As in 
most habitat types, vegetation plays a major role 
in determining the physical parameters of the 
habitat, and therefore, has a significant influence 
on the distributions of invertebrate species (Bonte 
2002, Tews et al. 2004). Gravesen (2000) and 
Lambeets et al. (2006) also support the hypothe- 
sis that a combination of flooding intensity and 
vegetation structure is important factor for the 
species composition of spiders in wet grassland 
areas. 

Sabo et al. (2005) reviewed the literature on 
the species richness and assemblage composition 
of flooded and non-flooded habitats. In accor- 
dance with the results of the present study they 
found that the species composition differs be- 
tween the regularly flooded and non-flooded hab- 
itats. This can increase the regional species rich- 
ness. Several studies have confirmed the impor- 
tance of fluvial dynamics affecting spider assem- 
blages by altering the habitat conditions (Bonn et 
al. 2002, Paetzold et al. 2008). 

At regularly flooded habitats Lambeets et al. 
(2009) found that the dispersal ability of species 
plays a prominent role in structuring carabid and 
spider assemblages, thus life-history traits affect 
species composition of spiders. Rothenbticher 
and Schaefer (2006) studied the impact of flood- 
ing on invertebrate assemblages. They distin- 
guished two types of adaptation in invertebrates: 
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submersion tolerance and regular migration be- 
fore and after the flooding. Most spiders over- 
winter at juvenile or adult stages (Pekar 1999). 
These stages of terrestrial arthropods are less 
flood-resistant than are their eggs, thus they are 
not able to tolerate long winter submersions. 
Rothenbiicher and Schaefer (2006) demonstrated 
that the floodplain spider fauna consists of spe- 
cies immigrating with receding water levels. The 
abundance of cursorial spiders depends on the 
proximity and connectivity of natural habitats, 
possibly for seasonal migration towards hiberna- 
tion sites (Bonte et al. 2003, Lambeets et al. 
2008), while good ballooners are able to colonize 
form greater distances (Bonte et al. 2002). 

Although the IndVal analysis identified 14 
species associated with floodplain grasslands, 
many species of the floodplain habitats are oppor- 
tunistic species from the surrounding non- 
flooded habitats that can colonize after the flood- 
ing events (Adis & Junk 2002). 


4.3. The proportion of forests 
had a significant negative effect 
on the rarefied number of the species 


In accordance with the results of the present 
study, Clough et al. (2005) reported that the rar- 
efied number of species is influenced by the land- 
scape parameters. Oberg et al. (2007) and Drape- 
la et al. (2008) found that a higher proportion of 
forests in the surrounding landscape were associ- 
ated with a higher number of observed spider spe- 
cies, as forests may serve as source habitat for spi- 
ders. Our results also support this hypothesis as 
we found increasing number of forest specialist 
species with the proportion of surrounding fo- 
rests. Species arriving in grasslands from forests 
are likely to add new species to the spider assem- 
blage, as forests and open habitats are home to 
contrasting spider assemblages (Entling et al. 
2007, Drapela et al. 2008). 

The effect of forests on the grassland special- 
ist species is possibly due to the fact that smaller, 
more isolated fragments retain fewer habitat spe- 
cies than larger less isolated ones (e.g. Lövei 
2006). Several papers emphasize the effects of 
surrounding habitats on the dispersal of species 
(Bonte et al. 2004, Thorbek & Topping 2005) and 


thus the isolation of the habitats. Higher propor- 
tion of forests may serve as dispersal barriers for 
grassland specialist species. 


4.4. The fauna of the surrounding forests 
influence the species composition of spiders 


According to the MRT landscape composition is 
important in determining the species composition 
in grasslands. Similarly to our results, numerous 
studies have shown that landscape variables are 
important determinants of species composition of 
spiders (e.g. Oberg et al. 2007, Batary et al. 2008, 
Schmidt et al. 2008). The IndVal analyses identi- 
fied two forest specialist species (P. alacris and 
D. concolor) at non-flooded grasslands with a 
high proportion of surrounding forests, and ac- 
cording to the linear models the number of forest 
specialist species increased and the number of 
grassland specialist species decreased with the 
proportion of surrounding forests. Although 
Schmidt et al. (2005) found no significant corre- 
lation between the percentage of forests and spi- 
der assemblages, Kajak (2007) demonstrated that 
forests influence the species composition of ce- 
real fields by enhancing the number and abun- 
dance of forest specialist species. She also found 
significant effect of forest age and naturalness. 

Landscape scale forestry operations may af- 
fect the diversity and species composition of the 
surrounding grassland habitats by altering the 
proportion or the habitat quality of forests. This 
emphasizes the importance of the landscape-level 
approach to nature conservation. 


5. Conclusions 


In the present study we have assessed the compo- 
sitional differences of spider assemblages of vari- 
ous grassland types and the effects of habitat and 
landscape properties on the rarefied species rich- 
ness and composition of spiders. Lowland flood- 
plain habitats are of high conservational value, as 
their associated arthropod fauna consists of spe- 
cialist and rare species (Rothenbiicher & Schae- 
fer 2004). Even small fragments play an impor- 
tant role in maintaining the regional diversity of 
spiders, as they are not entirely isolated, because 
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spiders are especially successful in dispersal 
(Horvath et al. 2009). 

A total number of 95 species were collected, 
and the species richness of grassland patches 
ranged between 23 and 34, thus the species com- 
position differed considerably between the stud- 
ied grasslands. Our results suggest that beside the 
conservation of floodplain grasslands, the main- 
tenance of diversity of grassland habitats is also 
essential for preserving the regional species pool 
of invertebrates. 
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